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Agitated  sand  bed 

Bottom  friction 

Energy  dissipation 
Nearshore  waves 

Nonlinear  wave  propagation 

Rough  turbulent  flow 

Have  height 

Wave  shoaling 

— >  This  report  provides  e  simplified  calculetlon  procedure  for  nearshore 
wave  height  changes  considering  the  energy  dissipated  by  rough  turbulent  flow 
over  a  strongly  agitated  bed  of  quarts  sand.  All  alien tary  wave  relation¬ 
ships  are  fron  linear  nonochronatlc  wave  theory,  but  one  effect  of  Including 
dissipation  is  that  calculated  height  changes  depend  on  the  absolute  wave 
height.  The  general  effect  of  appreciable  energy  loss  is  to  sake  field  wave 
height  relatively  constant  outside  the  breaker  sons.  Example  confutations 
and  a  calculator  progran  are  provided . "*r~  _  _ 
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this  report  provides  s  calculation  procedure  for  nearshore  shoaling  of 
energetic  neves  outside  the  breaker  cone,  including  the  appreciable  effects  of 
energy  dissipation  over  e  strongly  agitated  send  bed.  The  work  reported  was 
conducted  under  the  D.S.  Any  Coastal  Engineering  Research  Center's  (CERC) 
Numerical  Modeling  of  Shoreline  Response  to  Coastal  Structures  work  unit.  Shore 
Protection  and  Restoration  Research  Program,  Coastal  Engineering  Area  of  Civil 
Works  Research  and  Development. 

The  present  treatment  replaces  guidance  previously  provided  in  CERC  Field 
Guidance  Letter  Mo.  79-04  and  CERC  Technical  Paper  Mo.  80-S  (Grosskopf,  1980). 
Those  publications  incorrectly  recnund  the  use  of  friction  coefficients  for 
inert  beds  in  computing  nearshore  wave  shoaling. 

The  report  was  written  by  Dr.  Robert  J.  Hall erne lar.  Oceanographer ,  under 
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common  FACTOBS,  U.S.  CUSTOKABX  TO  WTBIC  (SI)  OUTS  OF  KASUBBOH 


U.S.  cuatomary  units 
aetric  (SI)  units  as 

of  aeasureaent 
follows: 

used  in 

this  report  can  be  converted  to 

Multiply 

ky 

To  obtain 

inches 

25.4 

milllmetera 

✓ 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

maters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

■lies 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

IQ'3 

kilograms  per  square  centimeter 

ounces 

28.35 

grama 

pounds 

453.6 

grama 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

1Tb  obtain  Celsius  (C)  temperature  readings  fro*  Fahrenheit  (F)  readings, 
use  formula:  C  »  (5/9)  (F  >32). 

To  obtain  Kelvin  (K)  readings,  use  foraule:  K  -  (5/9)  (F  -32)  ♦  273.15. 
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e  *  L/T  wr«  celerity 

D  asdlan  send  grain  disaster 

d  seen  water  depth 

dg  nartwua  water  depth  for  send  bed  agitation  by  waves 

K  average  energy  dissipation  rate 

f,  energy  dissipation  coefficient  for  rough  turbulent  flow  over 

strongly  agitated  sand  bed 

g  acceleration  due  to  gravity 

H  wave  height 

equivalent  wave  height  In  deep  water  Ignoring  refraction 
Ks  shoaling  coefficient  In  linear  wave  theory 

L  wavelength 

L,,  -  gT2/2*  wavelength  In  deep  water 
n  ratio  of  group  velocity  to  wave  celerity 

P  average  wave  energy  flux 

T  wave  period 

X  wave  propagation  distance 

€  horizontal  amplitude  of  near-bed  fluid  orbit 


p  fluid  density 

Additional  subscripts 

J  value  at  location  where  wave  condition  Is  to  be  predicted 

•  value  at  geoaetrlc  aeen  weter  depth  for  region  of  Interest 

1  value  at  location  where  wave  condition  Is  defined  initially 

+/-  case  where  dj  Is  greeter  or  less  than  dj 


CALCULATION  OF  WAVE  SHOALING  WITH  DISSIPATION 
OVER  NEARSHORE  SANDS 


by 

Robert  J.  Ballemeier 
I.  INTRODUCTION 

As  waves  propagate  toward  breaking  In  shallow  water,  their  attributes  are 
transformed  by  effects  of  water  depth  and  bottom  features.  These  nearshore 
transformations  are  crucial  in  the  interpretation  of  wave  measurements  and  in 
the  prediction  of  sediment  transport.  Prior  to  wave  breaking,  appreciable  wave 
energy  can  be  dissipated  by  friction  between  the  oscillatory  water  motion  and 
the  nearshore  bottom,  especially  where  waves  cause  strong  agitation  of  bottom 
sediments.  This  report  considers  such  frictional  dissipation. 

Ocean  waves  can  be  represented  most  adequately  as  distributions  of  propa¬ 
gating  energy  with  respect  to  frequency  and  direction.  Near  the  shore, 
extremely  energetic  waves  are  observed  to  constitute  a  somewhat  regular  wave 
train  approaching  along  a  shore-normal  line,  with  a  rather  well-defined  wave 
height  and  period.  According  to  the  Shore  Protection  Manual  (SIM)  (U.S.  Any, 
Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977),  an  idealisation 
of  demonstrated  value  in  coastal  engineering  is  to  represent  reel  waves  by  a 
characteristic  height  and  period,  the  significant  wave  condition;  this  wave 
representation  is  utilized  here. 

This  report  provides  a  simple  calculation  procedure  defining  changes  in 
significant  wave  height  due  to  water  depth  differences  combined  with  bottom 
friction  at  agitated  sand  beds.  The  procedure  uses  linear  wave  theory  between 
the  two  depths  of  interest  and  a  computation  of  energy  dissipation  rate  at  an 
Intermediate  water  depth.  Factors  in  wave  transformation  ignored  here  Include: 
wave  direction  and  complex  spectra,  currents,  surface  wave  breaking,  winds, 
water  viscosity,  and  bottom  percolation  and  elasticity.  The  technique  is  meant 
for  application  only  to  energetic  field  wave  conditions  in  fslrly  shallow  water 
with  straight,  parallel  depth  contours  and  relatively  fine  quartz  bottom  sands. 

Section  II  presents  the  calculation  procedure  for  wave  height  changes 
considering  energy  dissipation.  Section  III  addresses  the  application  of  this 
procedure  and  includes  two  example  problems:  converting  nearshore  wave  measure¬ 
ments  into  equivalent  wave  heights  in  shallower  and  in  deeper  water.  The 
reader  is  referred  to  Hallermeler  (in  preparation,  1983)  for  a  detailed  sub¬ 
stantiation  of  elements  incorporated  in  the  calculation  procedure;  that  refer¬ 
ence  reviews  the  empirical  basis  for  the  expression  giving  energy  dissipation 
coefficient  and  reports  extensive  calculated  results  in  clear  agreement  with 
multiple  wave  measurements  at  the  Coastal  Engineering  Research  Center's  (CERC) 
Field  Research  Facility,  Duck,  Nor*^  Carolina,  and  at  other  sites. 

II.  CALCULATION  PROCEDURE 

Waves  are  presumed  to  travel  perpendicular  to  depth  contours  with  propaga¬ 
tion  described  by  small-amplitude  (linear)  wave  theory.  Required  theoretical 
relations  are  provided  in  Section  2.23  of  the  SFM.  Wave  and  fluid  character¬ 
istics  arising  in  the  wave  description  are 
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d  «  local 


water  depth 


g  •  acceleration  due  to  gravity 
H  ■  local  wave  height 
L  ■  local  wavelength 
T  •  wave  period 
p  “  fluid  density 

Wavelength  in  deep  water  is  L©  ■  (gT2/2w)  and  the  dimensionless  local  wave¬ 
length,  d/L,  la  the  solution  of 


which  is  presented  in  Table  C-l  of  the  SFM.  The  average  wave  energy  flux  per 
unit  crest  width  is 

P  ■  |  P  g  H2  c  n  (2) 

where  c  -  (L/T)  la  wave  celerity  and  n  the  ratio  of  group  velocity  to  wave 
celerity: 

-  _  I  T,  ^  4rd/L  1 

(4ird/L)  J 


With  no  wave  refraction,  and  provided  no  energy  has  been  added  to  or  removed 
from  the  wave  train,  it  is  convenient  to  express  wave  height  changes  by  the 
factor 


where  la  equivalent  wave  height  in  deep  water  (ignoring  refraction),  and 
K8  the  shoaling  coefficient.  The  dimensionless  quantities  n  and  Kg  are 
provided  for  specific  values  of  (d/L0)  in  Table  C-l  of  the  SPM.  In  situations 
considered  here,  a  nearshore  wave  aeasursaient  at  water  depth  di  is  to  be 
converted  into  the  corresponding  wave  condition  at  another  relatively  shallow 
depth,  d j .  The  wave  period  is  presumed  constant  during  propagation  so  that 
(di/L0)  and  (dj/L0)  are  known,  and  equation  (4)  would  give  the  ratio  of  wave 
heights  as 

*.{*  <» 

“1  *S1 


if  energy  dissipation  were  to  be  ignored. 
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Energy  lost  from  the  ware  train  due  to  bottom  friction  is  treated  by  aeons 
of  a  single  dissipation  calculation  at  the  geometric  sen  depth  for  the  region 
of  Interest 


<*n  -  ^dj  <«> 

Average  energy  dissipation  rate  at  d*,  per  unit  crest  width  and  per  unit 
length  in  the  propagation  direction,  is  given  by 

-  0.235  p  feu  (2*  5u/T)3  (7) 

For  rough  turbulent  flow  over  a  strongly  agitated  bed  of  quarts  sand,  the  energy 
dissipation  coefficient  Introduced  in  equation  (7)  is 

fe»  -  exp  [-5.882  +  14.57  OWCn)°-m]  (8) 

Here  is  median  sand  grain  diameter  at  d«  and  Cm  1®  horizontal  ampli¬ 

tude  of  the  near-bed  fluid  excursion  arising  at  da  without  energy  dissipation, 
so  that  (2*€b/T)  in  equation  (7)  is  peak  near-bed  fluid  velocity.  According 
to  linear  wave  theory. 


Cm 


_ ^n _ 

2  slnh  (2rda/La) 


(9) 


where  Ha  “  (Hi  K3m/Kal )  from  equation  (5). 


Energy-conserving  linear  wave  shoaling  is  combined  with  computed  energy 
dissipation  rate  into  an  expression  giving  a  wave  height  at  dj  equivalent  to 
measured  wave  height  at  di.  This  expression  is  a  revised  form  of  equation 
(2): 


8(?i  ±  yo 

P  8  cj  “j 


(10) 


where  X  is  the  wave  propagation  distance  between  the  two  water  depths  (di 
and  dj),  and  the  upper  {lower]  sign  is  used  when  dj  is  greater  [less]  than 
dj.  The  conversion  given  in  equation  (10)  presumes  that  computed  dissipation 
rate  at  dm  can  be  considered  representative  of  the  entire  propagation  path. 


III.  APPLICATIONS 

Besides  the  explicitly  Ignored  factors  affecting  nearshore  wave  transfor¬ 
mations,  it  is  important  in  applications  to  consider  the  requirements  stated 
above  on  the  use  of  equation  (8)  for  energy  dissipation  coefficient.  The 
quarts  sand  bed  must  be  strongly  agitated  by  wave  action  and  near-bed  flow 
must  be  rough  turbulent.  Appropriate  situations  correspond  to  nearshore  field 
waves  with  relatively  large  height  and  period. 


The  strength  of  bed  agitation  nay  be  judged  using  an  approximate  expression 
(Hallemeiar,  1981,  eq.  10)  giving  nsxlawm  mater  depth,  d*.  for  wave  agitation 


(11) 


of  a  quarts  sand  bad  whan  viscous  offsets  are  negligible: 

d*  -  HT  (g/5000  D)°‘s 

If  da  computed  from  Hl(  Tlt  and  By  Is  much  larger  than  the  as*  tan  water 
depth  of  Interest,  the  requlreaent  for  strong  bed  agitation  nay  be  considered 
satisfied. 

To  assess  whether  flow  Is  likely  to  be  rough  turbulent,  another  simple 
computation  can  be  performed.  Incorporating  the  saae  approximation  for  £  in 
intermediate  water  depth  (2*d/L  near  unity)  as  is  utilized  in  equation  (11), 
fundamental  results  reviewed  in  Hallermeler  (in  preparation,  1983)  support 

HT  _>  d  {metric  units}  (12) 

as  an  approximate  criterion  for  rough  turbulent  flow  at  a  strongly  agitated  bed 
of  quarts  sand.  If  equation  (12)  Is  true  according  to  H},  Tj,  and  the 
!«<««■  water  depth  of  interest,  the  requirement  for  rough  turbulent  near-bed 
flow  nay  be  considered  satisfied. 

The  following  example  problems  demonstrate  the  use  of  the  present  procedure 
in  calculating  nearshore  wave  shoaling  with  energy  dissipation  due  to  a  strongly 
agitated  sand  '?ed . 


***************  example  PROBLEM  1*************** 

GIVEN:  At  the  CERC  Field  Research  Facility,  significant  wave  height  exceeding 
3.5  meters  was  recorded  during  three  1981  storms  by  a  Waverlder  buoy  located 
in  an  18-meter  Man  water  depth.  Have  periods  associated  with  these  ex¬ 
tremely  high  waves  ranged  from  9.3  to  14.0  seconds.  Nearshore  bathymetry  at 
this  site  is  regularly  surveyed  to  the  9-meter  water  depth  contour,  and  the 
wave  characteristics  at  the  seaward  Boundary  to  the  survey  region  are  of 
Interest.  The  shore-normal  distance  between  the  water  depths  is  1800  meters, 
and  the  representative  sand  slse  for  the  intervening  bottom  is  D  ■  0.12 
millimeter. 

FIND:  T.w  wave  height  at  dj  •  9  meters  corresponding  to  H|  -  3.5  Mters  at 
d]  -  18  meters  for:  (a)  T]  -9.3  seconds  and  (b)  Tj  -  14.0  seconds. 

SOLUTION: 


(a)  For  Tj  ■  9.3  seconds  and  d}  -  18  meters 

!i  .  "  (2w)(l»  -  0.1333 

lo  "  5tT  9.81  (9.3)2 

Table  C-l  in  the  SPM  gives  nx  -  0.7570,  (Hj/H*)  -  0.9160  -  It,}  and 

di/Lj  ■  0.1694,  so  that  Li  ■  106.3  maters.  With  Hi  ■  3.5  meters  and 
p  ■  1026  kilograms  per  cubic  meter  for  saltwater,  equation  (2)  becomes 


m 


fl  -  |  P  g  Hf  ci  nj  -  |  (1026) (9.81) (3. 5) 2  ^§2.  C0.757) 


-  1.33  *  10s  kilogram-meter  per  second  cubed 
Since  (Iq  my{d^  -  ^  (9) (18)  -  12.73  neters 


^m  (2w) (12.73) 
t  "  9.81  (9.3)2 
o 


0.09427 


and  Table  C-l  gives  (djj/I^,)  “  0.1360,  sinh  (Iwdm/l^)  “  0.9621,  (Hgi/Ho)  - 
0.9378  -  Ksm,  and  a*  -  0.8199.  Thus,  H»  -  (H,  1,,/!,^  -  (3.5) (0.9378)/ 
(0.9160)-  3.58  meters  and  according  to  equation  (9), 


With  Dm 


(m 


270621) 


1.86  meters 


0.12  millimeter,  equation  (8)  Is 


fern  -  exp  [-5.882  +  14.57  (Ve.) °* m] 

-  exp  [-5.882  +  14.57  (1.2  •  10'4/1.86)°*19,f]  -  0.0262 
so  that  equation  (7)  becomes 

14 . 0.235  p  “  (0.235)  (1026)  (0.0262)  [2,r|i^6lj3 

»  12.5  kilograms  per  second  cubed 
Some  conditions  must  also  be  computed  at  dj  —  9  meters,  where 


(2w)  9  « 

— J-  m  —  -  0.06665 

Lq  (9. 81) (9. 3) 2 

so  that  Table  C-l  gives  (H./H^)  -  0.9779  -  Kgj,  nj  «  0.8688,  and  (dj/Lj)  - 
0.1108  so  that  Lj  -  81.2  meters.  Finally,  because  dj  <  dj  the  lower  sign 
In  equation  (10)  is  appropriate,  and  X  -  1800  meters  yields 

„2  .  8(?l  -  V>  .  »H-33  •  105  -  (12.5)  (HOP)  |  .  u  5>  ,qu,t, 

3  P  8  CJ  (1026) (9.81)  g /3~  (0.8688) 

Hj  -  3.40  meters 

Prom  equation  (11),  maximum  water  depth  for  bed  agitation  is  da  ■ 
HjTj(g/5000D)®*5  -  (3. 5)(9. 3)[9.81/(5  •  103)(0.12  •  10“3)]0*5  -  131.6  meters. 


It 


much  larger  than  water  depths  In  the  region  treated,  and  the  numerical  value 
in  metric  units  of  (H jT j)  -  32.6,  nearly  twice  the  maximum  water  depth  con¬ 
sidered  in  meters,  so  that  equation  (12)  Indicates  rough  turbulent  flow 
throughout  the  region.  The  calculation  procedure  is  suitable  for  these  con¬ 
ditions,  and  the  effect  of  bottom  friction  on  wave  shoaling  is  appreciable, 
in  that  linear  wave  theory  without  dissipation  would  predict  a  nearshore  wave 
height  of  (Hi  Kgj/Ksl)  ■  [(3.5) (0.9779)/0. 9160]  ■  3.74  meters,  using  equation 
(5) . 

(b)  For  Ti  “  14.0  seconds  and  dj  ■  18  meters,  (di/L.)  *  0.05882  so  that 
«1  -  0.8833,  (Hj/H’)  -  0.9963  -  Ksl,  (dj/L^  -  0.1031  and  Lx  -  174.6  meters 

from  Table  C-l.  Equation  (2)  with  Hi  *  3.5  meters  gives 


Pj  -  |  (1026)  (9.81)  (3. 5)  2  (0.8833) 

■  1.70  •  10 5  kilogram-meter  per  second  cubed 

At  d_,  (djn/Lo)  -  0.04160  so  that  (d^L,,)  -  0.08509,  sinh  (2irdm/Lm)  -  0.5605, 
(Hn/Hj)  -  1.057  -  K^,  and  n  -  0.9161.  Thus,  H_  -  (H,  K^/Km)  -  3.71  meters, 
equation  (9)  is 


3.71 

2(0.5605) 


3.31  meters 


equation  (8)  is 

fem  "  e*P  (-5.882  +  14.57  (1.2  •  10-4/3.31)°* 194}  -  0.0207 
and  equation  (7)  is 

\  -  (0.235)  (1026)  (0.0207)  ]3 

*  16.36  kilograms  per  second  cubed 

At  dj,  (dj/L*,)  -  0.02941  so  that  (Hj/H<5)  -  1.130  -  Ksj,  nj  -  0.9400, 

Wj/Lj)  ■  0.0706  and  Lj  ■  127.5  meters.  Again,  the  lower  sign  in  equation 
(10)  is  used  to  give 

8[1.70  •  105  -  (16.36) (1800)] 

— * - * - - u—  -  13.05  square  meters 

(1026)  (9.81)  (0.9400) 

Hj  -  3.61  meters 

Because  T  is  greater  here  than  in  part  (a),  the  requirements  for  rough 
turbulent  flow  at  a  strongly  agitated  sand  bed  are  clearly  satisfied.  Although 
the  computation  including  friction  results  in  Hj  slightly  larger  than  Hj 
in  this  case,  energy  dissipation  again  has  an  appreciable  effect  since  linear 
wave  theory  would  predict  a  nearshore  wave  height  of  (Hj  K8j /K8j)  “  3.97  meters, 
from  equation  (5) . 
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***************  example  PROBLEM  2*************** 

GIVEN:  A  mathematical  model  Is  to  be  used  to  simulate  storm  wave  effects  for 
water  depths  shoreward  of  9  meters  at  Nags  Head,  North  Carolina,  with  the 
threshold  for  storm  waves  taken  to  be  the  wave  height  exceeded  10  percent  of 
the  time.  The  wave  climate  at  this  site  has  been  defined  by  several  rela¬ 
tively  complete  years  of  data  from  a  pier-mounted  gage  located  In  mean  water 
depth  of  5.2  meters  (Thompson,  1977):  wave  height  exceeded  in  10  percent  of 
these  measurements  is  about  1.7  meters  and  the  typical  wave  period  for  this 
wave  height  is  about  8.5  seconds.  The  shore-normal  distance  between  5.2-  and 
9-meter  water  depth  is  600  meters;  representative  sand  size  for  the  inter¬ 
vening  bottom  is  D  ■  0.20  millimeter. 

FIND:  The  wave  height  at  dj  -  9  meters  corresponding  to  Hj  «  1.7  meters  and 


Ti  “  8.5  seconds  at  d 


t  di 

1  ”  5. 


2  meters. 


SOLUTION:  For  Tj  -  8.5  seconds  and  d}  ■  5.2  meters. 


dj  2wd1  (2ir)  (5 . 2) 

L^  ’  gT Y  "  (9.81) (8.5)2 


0.04610 


Table  C-l  in  SPM  gives  nj  -  0.9074,  (Hj/H^)  -  1.038  -  KS1  and  dx/Li  - 
0.09002,  so  that  Li  *  57.76  meters.  With  Hi  *  1.7  meters,  equation  (2)  is 

Pj-ipgHl  Cjnj-I  (1026)  (9. 81)  (1.7)  2  (0.9074) 

■  2.24  •  104  kilogram-meter  per  second  cubed 

Since  dm  -  /djdj  -  /(5.2) (9)  -  6.84  meters 

f«  .  (2*) (6.84)  m 

Lq  (9.81) (8.5) 2  °-0606A 

and  Table  C-l  gives  (dm/Lo)  -  0.1049,  sinh  (2irdm/l«)  -  0.7082,  0WH£)  - 
0.9916  -  Kg_,  and  n,,  -  0.8799.  Thus,  Hj,  -  (^  K^/K^)  -  (1.7)  CO. 9916)/ 

1.038  «  1.624  meters  and  according  to  equation  (9), 


2  sinh 


(*?) 


1.624 

2(0.7082) 


1.15  meters 


With  Dm  ■  0.20  millimeter,  equation  C8)  is 

fem  -  «P  (-5.882  +  14.57  (Dm/C*)0-194} 

-  exp  (-5.882  +  14.57  (2  •  10"4/1.15]°-194}  -  0.0422 


so  that  aquation  (7)  la 


\  -  0.235  p  tjfif*)  3  *  CO . 235) (1026) (0 .0422)  [—  fe1-5*-]* 

-  6.25  kilograms  pet  second  cubed 
At  dj  -  9  meters, 

il  -  — - -  0.0798 

^  (9.81) (8.5) 2 


so  that  Table  C-l  gives  (H«/H£)  ■  0.9551  »  Kg* ,  nj  »  0.845,  and  (dj/Lj)  * 
0.1230  so  that  Lj  •  73.2  asters.  Finally,  because  dj  >  dj  the  upper 
sign  In  equation  (10)  is  appropriate,  and  X  -  600  aeters  yields 


„2  .  8  (P»*M)  .  8  :  10“  *  <M5)WP>1  .  2.66  square  -ter. 

1  0  *  CJ  "J  (10 26) (9. 81)  i&Q  (0.8*5) 


(0.845) 


Hj  -  1.69  meters 

To  show  that  the  calculation  procedure  is  suitable  for  these  conditions, 
maximum  water  depth  for  bed  agitation  from  equation  (11)  is 


(1.7) (8. 5) 


[~ 


103) (2  •  10“4) 


->45.3  meters 


much  larger  than  water  depths  In  the  region  treated,  and  the  numerical  value  in 
metric  units  of  (HjT^  -  14.45,  greater  than  the  maximum  water  depth  considered 
in  meters,  so  that  equation  (12)  Indicates  rough  turbulent  flow  throughout  the 
region  considered.  The  effect  of  bottom  friction  is  still  appreciable  for 
this  relatively  low-energy  case,  in  that  linear  wave  theory  without  dissipation 
would  provide  a  wave  height  according  to  equation  (5)  of  (Hj  Kgj/Kg^)  ■ 
f(1.7)(0.9551)/1.038]  -  1.56  meters  st  9-meter  water  depth. 


***************  ************»»«*»,**»',',‘■c,* 

With  linear  wave  theory,  the  height  change  between  two  water  depths  depends 
on  wave  period.  Although  only  linear  theory  wave  relationships  are  Incorpora¬ 
ted  in  the  present  calculation  procedure,  energy  dissipation  depends  both  on 
wave  period  and  wave  height  (raised  to  the  power  of  about  2.5).  Thus,  the 
calculated  results  have  a  nonlinear  dependence  on  wave  height:  the  computed 
height  change  between  two  water  depths  Is  affected  by  the  actual  value  of  wave 
height. 
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This  non  Ilnur  aspect  lapliu  that  than*  computations  ara  not  aactly 
reversible .  Projecting  •  nave  condition  offshore  without  dissipation  from  a 
measurement  sits  to  ^  can  result  in  a  markedly  different  conputed  dissipa¬ 
tion  rate  there  than  if  the  nominally  corresponding  waves  are  projected 
onshore  to  dB.  However,  the  calculation  procedure  tends  to  cancel  Internally 
this  effect  of  nonlinear  height  dependence. 

Sever sing  Example  Problea  2,  using  Hj  -  1.69  asters  at  dj  •  9  asters  as 
the  Initial  condition,  coaputed  conditions  at  dm  include  £■  -  1.24  asters, 
fea  *  0.0406  and  Ea  -  7. 52  kilograms  per  second  cubed,  but  at  dj  ■  5.2  asters 
the  calculated  wave  height  is  1.67  asters,  only  1.8  percent  less  than  the  near¬ 
shore  wave  height  df  1.7  aetera  originally  specified.  Using  coaputed  final 
wave  heights  in  Example  Problem  1  as  input  conditions,  the  reverse  calculation 
procedure  gives  wave  heights  in  each  case  only  .2. percent  less  than  the  specified 
height  of  3.5  meters. 

Such  slightly  irreversible  results  do  not  seem  too  significant  for  potential 
applications.  However,  the  Appendix  to  thia  report  provides  a  calculator 
prograa  quickly  executing  the  present  procedure,  asking  it  convenient  to 
examine  results  of  reverse  calculations  and  to  determine  a  wave  condition  which 
appears  optimally  consistent  with  that  specified. 

IV.  SUMMARY 

The  equations  and  procedures  presented  here  permit  calculation  of  nearshore 
wave  height  changes  considering  the  energy  dissipated  by  rough  turbulent  flow 
over  a  strongly  agitated  bed  of  quarts  sand.  All  elementary  wave  relationships 
are  from  linear  (small-amplitude)  wave  theory,  but  one  effect  of  incorporating 
dissipation  is  that  calculated  height  changes  depend  on  the  actual  wave  height. 
Example  calculations  demonstrate  the  conversion  of  a  nearshore  wave  condition 
into  a  corresponding  wave  height  in  shallower  or  deeper  water;  the  present 
procedures  are  suitable  only  for  field  waves  of  relatively  large  height  and 
period  In  fairly  shallow  water.  The  general  effect  of  energy  dissipation  ia 
that  nearshore  wave  height  remains  more  nearly  constant  outside  the  breaker 
sone  than  linear  wave  theory  would  predict. 
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APPENDIX 


CALCULATOR  PROGRAM  FOR  WAVE  SHOALING  WITH  DISSIPATION 

The  following  four  pages  docuaont  a  calculator  prograa  executing  the 
procedure  presented  in  Section  II  and  exemplified  in  Section  III  of  this  report. 
This  prograa  runs  In  about  120  seconds  on  a  Hewlett-Packard  HP-67  Programmable 
Pocket  Calculator,  employing  aetrlc  units,  KPN  logic,  three  levels  of  nested 
subroutines,  17  address  labels,  26  storage  registers,  and  223  prograa  steps. 

The  prograa  could  be  converted  for  use  with  other  calculators  having  different 
logic  systems  but  similar  features  and  capacities. 

Equations  (1)  to  (10)  are  Included  with  an  effective  root-finding  Iteration 
for  wavelength.  Values  to  be  specified  In  metric  units  for  a  calculation  are: 
p,  g,  Hj,  Tj,  dj,  dj,  D#,  and  X.  The  standard  value  of  g  Is  9.81  maters  par 
second  squared,  and  the  value  of  p  for  seawater  nay  be  taken  as  1026  kllo- 
graas  per  cubic  aeter;  for  freshwater,  p  la  about  1000  kilograms  per  cubic 
aeter,  but  coamon  situations  night  not  constitute  the  requisite  rough  turbu¬ 
lent  flow  over  a  strongly  agitated  bed.  Satisfaction  of  these  requirements, 
related  to  equations  (11)  and  (12) ,  remains  to  be  considered  external  to  the 
calculator  prograa. 
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